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ABSTRACT 



o 



. We determine the metallicity distribution function (MDF) of the Galactic halo based on metal-poor main-sequence turnoff-stars 
Q_l' (MSTO) which were selected from the Hamburg/ESO objective-prism survey (HES) database. Corresponding follow-up moderate- 
resolution observations (R - 2000) of some 682 stars (among which 617 were accepted program stars) were carried out with the 2.3m 
telescope at the Siding Spring Observatory (SSO). Corrections for the survey volume covered by the sample stars were quantitatively 
estimated and applied to the observed MDF. The corrections are quite small, when compared with those for a previously studied 
sample of metal-poor giants. The corrected observational MDF of the turnoff sample was then compared with that of the giants, as well 
as with a number of theoretical predictions of Galactic chemical evolution, including the mass-loss modified Simple Model. Although 
the survey-volume corrected MDFs of the metal-poor turnoff and the halo giants notably differ in the region of [Fe/H] > -2.0, 
, below [Fe/H] ~ -2.0, (the region we scientifically focus on most) both MDFs show a sharp drop at [Fe/H] ~ -3.6 and present 
■ rather similar distributions in the low-metallicity tail. Theoretical models can fit some parts of the observed MDF, but none is found 
to simultaneously reproduce the peak as well as the features in the metal-poor region with [Fe/H] between -2.0 to -3.6. Among 
the tested models only the GAMETE model, when normalized to the tail of the observed MDF below [Fe/H] ~ -3.0, and with 
. Z„. = 10~ 34 Z o , is able to predict the sharp drop at [Fe/H] 3.6. 
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1 . Introduction the early stages of galaxy evolution. Although we are gaining a 

deeper understanding over time, much remains to be explored. 
The Galactic halo provides important clues for understand- It is particularly revealing that, after many decades of assum- 
ing the evolution and structure of the Galaxy. In the past few ing that the Galactic halo compr ises a single stellar population, 
decades, considerable observational and theoretical efforts have recent work dCarollo et al.ll2007l (2010) has provided additional 
H , been made to investigate its chemical evolution, details of its support to suspicions mat emerged from previous efforts that the 
structure, and its kinematical characteristics. Very metal-poor halo ifj indeed divisible int0 two structura i components, with no- 
stars in the halo, those with metallicity [Fe/HlJ < -2.0, are re- tably differ ent spatial density profiles, stellar orbits, and stellar 
garded as fossils of the earliest generations of stars. They pre- metallicities 
serve the chemical information created by their stellar progeni- 
tors, providing fundamental insights regarding the properties of R ecently, new theoretica l models (e.g jHelmj2008b lP rantzosl 

the very first generation of stars, the chemical history of our 1 20081: ISalvadori et al1l2010l) and observational constraints (e.g. , 

Galaxy (and other large spirals like it), the modes of star for- I Carollo et all 120071 1201 Ot iBell et al. l j 2008; Ive zic et al.1 [2008; 

mation in the proto-Milky Way, the formation of the Galactic lJuric et alJ 120081: iBo nd et alJ 120091: Ide Jong et alJ I201QI) have 

halo, and physical mechanisms such as feedback processes in greatly enhanced our understanding of the nature of the halo 

components of our Galaxy, Those are enabling the development 

* Based on observations collected at Siding Spring Observatory. of plausible assembly histories based on the degree of detectable 

1 The common notation of [A/B] = log(N A /N B )*- log(N A /N B ) Q is spatial and phase-space substructures. The possible association 

used here, where N A and N B are the number densities of elements A and of at least some presently observed dwarf gala xies with the 

B, respectively. formation of the halo populations, as invoked bv ICarollo et al.l 
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(2007) to account for their dual halo structure, has received ad- 
ditional support based on high-resolution spectroscopic analy- 
sis of in dividual stars in ultra-faint and dwarf spheroidal galax- 
iesfe.g.. iMufioz et alJl2006l: iKirbv et alJl2008t iGeha et al.ll2009t 
iFrebel et alJl2010t lNorris et alj |20 1 Oh .Finally, the identification 
and detailed analysis of the elemental abundance patterns for the 
most chemically p rimitive stars, e.g., the ultra ( [Fe/Hl < -4.0; 
Norris et all 120071) and hyper ([Fe/Hl < - 5.0; IChristlieb et alj 
2002t IFrebel et al.ll2005t [Aoki et al.ll2006l) metal-poor stars al- 
low one to trace back close to the very beginning of star forma- 
tion in the Galaxy. 

The observed metallicity distribution function (MDF) of halo 
stars provides strong constraints on models for the formation 
and chemical evolution of the Galaxy. Any accepted model 
must be able to predict the r elative numbers of halo stars as 
a fun ction of their metallicity dBeers & Christliebl l200l iHermil 
2008), and in the case of a dual-halo model, as a function of 
location and kinematics. Early investigations on the shape of 
the halo MDF were hampered by the small numbers of very 
metal-poor stars known at the time dHartwickHl976t iBondf l 981; 



Rvan&Norrisl|199U 
Bonifacio et alkoOOt 



Carney et al. 1996). Other attempts (e.g., 
Schuster et al.ll2.004l) . based on samples of 
metal-poor stars from th e HK survey of Beers and colleagues 
dBeers et al.lll985l [l992). suffer from poorly constrained selec- 
tion criteria, except perhaps at the lowest metallicities. More re- 
cent efforts have made use of statistically well-understood selec- 
tion criteria to identify large numbers of metal-poor candidates 
from obje ctive-prism surveys , such as the Hamburg/ESO survey 
(HES -[Wi sotzki et al. 1996 ), as reported in a series of papers 
dBarklem et al Jl2005h IChr istlieb e taU20q8t ISchorck et alj|2009t 
Placc o et alJl2010h . ISchorck et alJ (l2009h . for example, used a 
sample of 1638 metal-poor giants to study the shape of the low- 
metallicity tail of the halo MDF, and made detailed comparisons 
with MDFs of Galactic globular clusters and satellite galaxies, 
as well as with theoretical models. 

Main-sequence turnoff (MSTO) stars have long been used 
to explore Galactic structure, including the recognition of 
stella r substructures in the Galactic halo dMaiewski et alj 
2004; Anetal. 2009), searches for kinematic streams (e.g., 
Klement et alJ l2009). and statistical analyses of the amount of 
cold halo substructure in the Milky Way (e.g. JSchlaufman et alj 
2009). In addition, MSTO stars have also been proven im- 
portant to the field of Galactic chemical evolution, through 
the analysis of high-resolution, high signal-to-noise spectro- 
scopic observation s to derive element al abundances for metal- 
poor dwarf stars dCohen et al. 2004)). ch emically interesting 
metal-poor turnoff stars dAoki et al.l 120081). and in vestigations 
of the so-called Spite Plateau dSpite & Spitdll982h through Li 
abundance measurements fo r metal -poor turnoff stars t Aoki et alJ 
2009; Sbo rdone et alj2010l) . In this paper we construct the MDF 
of Galactic halo MSTO stars based on follow-up moderate- 
resolution (R ~ 2000) spectroscopic observations of candidate 
metal -poor turnoff stars from the HES. We also compare our re- 
sults with the previously derived MDF of HES giants, and with 
theoretical expectations. 



This work is 
series (Pape r I - 



a continu ation of the 
Christlieb 



HES stellar content 
eb et alJ l200Tbl II - IChristlieb > et alJ 
l2001al III - ICh ristlieb et al.l I2005L IV - Ichristlieb et all 12008. 
V - ISchorck et al.ll2009l) 7We describe the selection of the HES 
turnoff sample in Section [2] with details of the metallicity de- 
termination and MDF construction given in Section [3] The ob- 
served HES MSTO MDF is compared with theoretical predic- 
tions in Section|4l the main results are summarized in Section[5] 



2. The Sample 

Adopting the methods described in Paper IV, the HES metal- 
poor turnoff candidates were selected from the HES objective- 
prism database, using both KP/(B - V)o and KP/(7 - K)q selec- 
tions. The only exception was that an additional (B - V)o range 
was specified so that the candidates were restricted to 0.3 < 
{B - V)q< 0.5. The KP index hereby m easures the streng th of 
the Ca ii K line and is defined in detail in Be ers et a D dl999h ). It 
was measured in all HES prism spectra and together with a color, 
the prime indicator for the selection of metal-poor candidates. 

As shown in Figure 6 of Paper IV, the employed KP cut- 
off becomes comparable to its measurement uncertainty for stars 
within our (B - V)o range. Since turnoff stars are also rela- 
tively weak-lined, we thus included only those candidates whose 
Can K line is not significantly detected in HES spectra, even 
if their KP indices are above the formal cutoff line. To main- 
tain relatively consistent exposure times during the follow-up 
observations, an additional cutoff regarding the brightness of 
Bhes ^ 16.5 was adopted. The above cuts yielded a prelimi- 
nary sample of 3383 metal -poor turnoff candidates from the HES 
database. 

In order to provide candidates with a higher likelihood of 
being metal-poor, the HES prism spectra of the selected 3383 
candidate were visually inspected. As defined in Paper IV, based 
on the apparent strength of the Ca n K line relative to the con- 
tinuum, the 3383 candidates were classified into four different 
metal-poor classes, mpca, unid, mpcb, and mpcc. The distribu- 
tions of the 3383 candidates for these classes are listed in the 
second column of Table [1] 

For accurate measurements of stellar metallicities (as well as 
the estimates of other stellar atmospheric parameters), moderate- 
resolution follow-up spectra are required. To avoid possible 
systematic offsets of spectral features that could arise from 
combining different telescope/detector combinations, we ex- 
clusively adopted data observed during 15 individual runs at 
the Siding Spring Observatory (SSO) 2.3m telescope with the 
Double Beam Spectrograph (DBS). The runs took place between 
January 2006 and November 2009. The resolving power was R 
- 2000, with a typical S/N of 20/1 per pixel in the continuum 
region close to the Can K line. 

For a total of 682 unique stars from our metal-poor turnoff 
candidate list follow-up spectra were obtained. The third column 
of Table Q] lists the numbers of these observed candidates for all 
four metal-poor classes. It is clear that a significant bias against 
the class mpcc exist because it is the subjectively least promising 
candidate class for finding metal-poor stars. 



Table 1. Numbers of different metal-poor classes among the can- 
didate HES turnoff sample, the sample with follow-up observa- 
tions from SSO, and the accepted sample after the visual inspec- 
tion and rejection as described in Section |3~T1 The column la- 
belled "Factor" refers to the scaling factor used to construct the 
MDF described in Sectionl3~Tl 



Class 


HES prism 


SSO follow-up 


Accepted 


Factor 


mpca 


179 


36 


29 


6.17 


unid 


333 


67 


59 


5.64 


mpcb 


1666 


560 


513 


3.25 


mpcc 


1205 


19 


16 


75.31 


Total 


3383 


682 


617 
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3. Analysis of the Observational Sample 

3.1. Measurements of line indices and [Fe/H] 

Line indices (Beer s et al. 1999) were measured for all 682 stars 
in our program sample which include the KP index, the HP2 in- 
dex (which measures the strength of the Balmer Hg line), and 
the GP index (which measures the strength of the of CH G- 
band feature). In the cases where we had multiple spectra for 
a single object, we adopted a S/N-weighted average of the in- 
dividual indices. Following this step, a visual inspection of the 
follow-up spectra was carried out to identify and reject spec- 
tra of objects that were too noisy, had emission lines present, 
or were too hot (as indicated by their Balmer line index); a few 
additional objects that turned out to be galaxies or were other- 
wise peculiar were rejected as well. No stars with GP > 6 A, 
which indicates strong spectral carbon features, were detected in 
our sample. This is perhaps not surprising, given the relatively 
high effective temperatures of our turnoff sample. We note that 
HE 1327-2326, the most metal-deficient star currently known, 
was one of the candidates observed during the 15 runs. On the 
grounds that the star was known to be hyper-metal-poor and in 
order to obtain a better medium-resolution spectrum than previ- 
ously existed, it was included in the follow-up observations de- 
scribed here. However, we removed it from our sample because 
including this star might introduce a bias to the sample. This left 
us with an "accepted" sample containing 617 stars, with the dis- 
tribution across different metal-poor classes given in Column 4 
of Table □ 

To obtain [Fe/H] estimates for the stars in our sample, 
two independent procedures were carried out. The first ob- 
tains [Fe/H] from the measured KP and HP2 indices by us- 
ing an up dated code version making use of the methods de- 
scribed by iBeers et alj (Q999) (which includes more calibration 
stars, and thus results in a better coverage of stellar param- 
eter space, especially in the lowest metallicity regime). This 
was the method used in Paper V when constructing the giant- 
star sample. The second method is a newly develop ed version 
of the S EGUE Stellar Param eter Pipeline (SSPP - iLee et alJ 
2008a1bl: IXilende Prieto et aT][2008h . The SSPP is the software 
tool used to obtain estimates of atmospheric parameters for stel- 
lar spec tra obtained during the course of the Slo an Digital Sky 
Survey dYork et all 2000; Ab azaiian et al.ll2009t) and its exten- 
sions, SDSS-II and SDSS-III. Recent experiments with spec- 
tral data with similar resolving power to SDSS spectra (R = 
2000) have indicated that the SSPP can provide useful esti- 
mates of parameters for non-SDSS data as well, as long as the 
wavelength coverage extends from roughly 3800 A to 5200 A. 
Slightly smaller wavelength ranges can still be used, but the ac- 
curacy of the derived parameters (in particular surface gravity) 
begins to suffer when the red limit is less than 5000 A, due to the 
loss of the Mg lb and MgH features, which provide enhanced 
sensitivity to estimates of log g. 

The non-SEGUE Stellar Parameter Pipeline (n-SSPP) takes 
as inputs user-supplied measurements or estimates of the 
Johnson V ma gnitude and B - V color, and/or a 2MASS 
dCutri et al.ll2.Q03h J magnitude and J - K color, all corrected for 
a user-specified level of absorption and reddening, along with 
a user-supplied estimate of the observed radial velocity. It then 
proceeds to determine estimates of the primary atmospheric pa- 
rameters (T e s, log g, [Fe/H]) and their estimated errors, as well 
as estimates of di stance, making u se of a subset of the proce- 
dures described in lLee et aTl d2008al) (i.e., those that can be made 
to work within the wavelength region covered by the input spec- 



trum). Note that it is not necessary that the input spectra be flux 
calibrated, nor continuum rectified. It is also not strictly nec- 
essary to supply input colors, since the n-SSPP makes internal 
estimates that can be used as needed, but due to possible degen- 
eracies in the derived parameters color information is certainly 
preferred. 

The n-SSPP was used to obtain atmospheric parameter esti- 
mates for our 617 accepted program stars. Input Johnson colors 
were taken from the estimates provided by the HES catalog (for 
V and B - V), or photome tric measurements for a small subset 
of the data available from IBeers e t al. (2007) or later (Beers et 
al., in prep.); J and J - K were taken from the 2MASS Point 
Sourc e Catalog, absorption corrected or de-reddened according 
to the lSchlegel et alJd 19981) dust maps. 



0.5 r 




-4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 

[Fe/H] Beeral999 

Fig. 1. Comparison of the metallicities ([Fe/H]) of the HES 
turnoff sample m easured with the updated version of the 
IBeers et al] (11999b method and the n-SSPP. The dashed line is 
the one-to-one line; typical errors for both measurements are 
shown in the upper-left corner. A slight offset between the meth- 
ods is clearly present. See text for discussion. 

The resulting estimates of [Fe/H] for the two methods are 
compared in Figure Q] Although th e two measuremen ts do not 
greatly differ, the typical error of the lBeers et al 1 (fl999l) determi- 
nation (0. 18 dex) is twice tha t of the n-SSPP est imate (0.09 dex). 
The determination based on Beer s et al.l (1 19991) exhibits an off- 
set in [Fe/H] compared with the n-SSPP of -0.19 ± 0.01, re- 
sulting in a metallicity distribution that reaches apparently lower 
[Fe/H]. In order to make the derived metallicity distribution as 
accurate as possible we have adopted the high-resolution mea- 
surements for the three candid ates that have been observed with 
high-resolution spectroscopy dCohen et al.l l2004h . As shown in 
Table [2] the metallicities obtained by the n-SSPP for these stars 
are closer to the values derived by the high-resolution analysis 
than are those from the Beers et al. approach. Considering the 
fact that the n-SSPP procedure delivers what is likely to be a 
more accurate estimate of [Fe/H], the following statistical dis- 
cussion will be based on this method. 

Distances to our sample stars were calculated by assigning 
the stars into various luminosity classes. This was done based 
on the surface gravity estimates derived by the n-SSPP. We con- 
sidered all stars with estimated log g > 3.5 to be likely dwarfs, 
those with 3.5 < log g < 3.0 to be turnoff stars, and those with log 
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Table 2. Comparisons of difference measurement techniques 
for [Fe/H]. Columns 2-4 resp ectively refer to v alues derived 
by the updated version of the iBeers et al.l ( 1999) method, the 
n-SSPP, and detailed analysis by Coh en et al.l d2004l) based on 
high-resolution spectroscopy. 



HES ID 


[Fe/H] | 


[Fe/H] 2 


[Fe/H], 


HE 0007-1832 


-3.32 


-2.98 


-2.65 


HE 0105-2202 


-2.87 


-2.86 


-2.55 


HE 1346-0427 


-3.32 


-3.57 


-3.40 



g < 3.0 to be subgiants and giants. The distribution of the sam- 
ple in the Z-R plane is shown in Figure [2] It indicates that our 
sample of turnoff candidates is indeed located within 2 - 3kpc 
of the Sun. This can be contrasted with a similar figure (Figure 
3) from Paper V which shows that the giants are located at much 
larger distances from the Sun. 




5 6 7 8 9 10 11 
R/kpc 



3.3. Construction of the MDF 

As discussed above, the MDF derived from our follow-up ob- 
servations contains a significant bias towards the more metal- 
deficient candidates, and must be taken into account to recover 
a reasonable representation of the "true" MDF. Therefore, we 
adopted the scaling factor procedure described in Paper V. For 
each metal-poor class, the MDF of the observed candidates is 
scaled by a factor calculated from the division of the total num- 
ber in the class by the observed number (as listed in the last 
column in Table [TJ. Then the scaled MDFs of the four classes 
are co-added to produce a general MDF for the entire HES can- 
didate sample. Similarly to Paper V, the main difference between 
the directly observed and the scaled MDF is the increasing ratio 
of the relatively metal-rich stars in the mpcb and mpcc classes. 
The normalized fraction of the scaled MDF is listed in the first 
column of Table [3] 

3.4. Selection fraction 

As pointed out in Paper IV and V, the combination of the KP 
index with (B - V)o or (J - K)o for the purpose to select 
metal-poor candidates in the HES has proven rather efficient. 
Following the metallicity distribution predicted by the Simple 
Model, we apply our quantitative selection criteria to a simulated 
sample of metal-poor stars. The results of the theoretical selec- 
tion fractions shown in Figure[4] The selection fractions for both 
(B - V)o and (J - K)o are shown. It is clear that the selection cri- 
teria are able to reject the majority of stars with [Fe/H] greater 
than -2.0. For both colors, a high completeness (up to almost 
100%) is reached for stars with [Fe/H] < -3.0. For (B - V) , the 
redder candidates exhibit a larger selection fraction (due to less 
contamination from hot stars among the bluer candidates). The 
selection fraction, however, does not differ much among the dif- 
ferent (J - K)q cutoffs. This is as expected since the blue cutoff 
in (J - K)q is already fairly red so that fewer hot candidates enter 
the sample. 



Fig. 2. Spatial distribution of the observed HES MSTO sample. 
R indicates the Galactocentric distance projected to the plane; Z 
indicates the distance above or below the Galactic plane. In this 
diagram, the Sun is assumed to be located at R = 8 kpc, Z = 0. 



3.2. Observational biases and selection effects 

Before comparing our observed MDF with other results or the- 
oretical predictions, it is necessary to address the biases and se- 
lection effects that are introduced through the survey itself or by 
our sample selection procedure. 

One notable selection bias comes from the metal-poor classi- 
fication. As shown in Figure[3]the follow-up observations clearly 
favor the best metal-poor candidates. The numbers in the panels 
indicate that the relative selection efficiency of extremely metal- 
poor (EMP) stars ([Fe/H] < -3.0) obviously decreases from a 
maximum of about 10% for the better candidate classes, mpca 
and unid, to the least likely class, mpcc. Note that the fraction 
of EMP stars in class mpca is somewhat lower than that in class 
unid, probably due to the fact that turnoff metal-poor stars are 
rather weak-lined, making the divisions between these classes 
rather difficult. The numbers of targets in mpcc is rather small 
(16), hence it is perhaps not surprising that no EMP stars were 
found in this class. 



3.5. Survey volume correction 

As pointed out in Paper V, for a magnitude-limited survey the 
relative survey volume explored by the observed stars differs 
with the stars' metallicities, which could also be readily inferred 
from Figure|5] Besides, as described in Section|2]and Table[TJ the 
HES follow-up procedure is basically a metallicity-biased sur- 
vey, which favors candidates with lower metallicities. Thus it is 
interesting to investigate to what extent this effect could impact 
our sample and the resulting derived MDF. Moreover, we aim at 
deriving a corrected MDF that is metallicity/volume- unbiased 
suitable for the comparison with other observational results and 
theoretical models. 

The basic idea of this correction is to derive the survey vol- 
ume for stars with different metallicities, referenced to a specific 
metallicity. Here we adopt [Fe/H] = -2.0, because it is near 
the peak [Fe/H] of our sample, and also close to the metallic- 
ity above which we expect the observed MDF to deviate from 
the "true" MDF due to metallicity selection bias. It is thus con- 
venient for later comparisons (the choice of a different refer- 
ence [Fe/H] will not strongly affect the relative fraction of 
each [Fe/H] bin of the corrected MDF). Based on the defini- 
tion of the survey volume, the corrected volume referenced to 
[Fe/H] = -2.0 in a specific (B - V)o bin can be directly esti- 
mated from V = iQ0-6(.My(ref)-My) As forthe tU rnoff sample, stars 
within a (B - V) and [Fe/H] bin could be either a MSTO star 
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Fig. 3. The observed MDFs of the different metal-poor classes. The numbers listed in each panel correspond to the number of EMP 
stars (the numerator) compared to the total number of stars in the corresponding class (the denominator). 
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1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



Total 

(B-V) = 0.3 
(B-V = 0.4 
(B-V) = 0.5 





Total 

(J-K) = 0.3 
(J-K) = 0.4 
(J-K) = 0.5 
(J-K) = 0.6 
(J-K) = 0.7 



-4.0 -3.5 



-3.0 -2.5 
[Fe/H] 



-2.0 -1.5 -1.0 -4.0 -3.5 



-3.0 -2.5 
[Fe/H] 



-2.0 -1.5 -1.0 



Fig. 4. Selection fractions of metal-poor candidates selected with the HES selection criteria as described in Paper IV. The two panels 
correspond to selection efficiencies using KP and (B - V)o (left) or (J - K)q (right). Different lines refer to different red cutoffs as 
shown in the legend; the solid lines refer to the total selection fractions. 



or a subgiant, which obviously explore different survey volumes. 
Therefore, another step in the correction is used to estimate the 
ratio of the MSTO stars to subgiants in the sample. Using the lu- 
minosity functions from the Y 2 isochrones and assuming an IMF 
slope of x = 1.35 (Salpeter index), for any specific [Fe/H] and 
(B - V)owe can obtain the number of stars per cubic parsec per 
absolute magnitude interval for both the MSTO and subgiant 
branches. Hence a relative density ratio of MSTO stars versus 
subgiants for the sample is obtained. Given the relative number 
of MSTO stars and subgiants in each [Fe/H] and (B - V)o bin, 
we can then obtain the corrected number of stars within a spe- 
cific [Fe/H] and (B - V)o bin by combining the volume and the 
fraction corresponding to the MSTO and subgiant stars. 



Based on this procedure, we derive the volume-corrected 
MDF of the sample and compare it with the observed one. This 
is shown in the left panel of Figure [6] As can be seen from in- 
spection of this figure, the survey-volume correction only very 
slightly affects the shape of the MDF. It mildly decreases the 
fraction of lower metallicity stars (referenced to [Fe/H] = -2.0) 
while slightly increasing the fraction at higher metallicity. This 
is not unexpected because our sample of turnoff stars occupy 
a relatively narrow range of (B - V)q near the blue end of 
the isochrones (see Figure [5j. Thus, their relative observational 
volumes for different metallicities or different branches on the 
isochrones (MSTO or subgiant) do not greatly differ. This is 
also supported by the spatial distribution of our sample shown 
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Fig. 5. Theoretical is ochrones from iDemarque etal . (2004), based on the Yonsei-Yale isochrones Version 2 dYi et alJ 120011: 
Demarqu e et al.l2004l) . with parameters shown in the left panel; the dash-dotted, solid, and dotted lines correspond to [Fe/H] = - 1 .0, 
-2.0, and -3.0, respectively. The color ranges of the HES turnoff sample are indicated by the vertical dashed lines. 



in Figure [2] The correction factors of each [Fe/H] bin for the 
MDF are listed in the third column of Table [3] and are also ap- 
plied to the corresponding [Fe/H] bins of the scaled MDF of the 
complete candidate sample of 3833 stars derived in Section [331 
(given in the last column of Table[3j. 

Table 3. The survey-volume correction for the observed HES 
turnoff MDF. The second column refers to the fraction of the 
scaled MDF; the third column refers to the factor arising from 
the survey-volume effect as described in Section l3~5l and the last 
column refers to the corrected fraction for each [Fe/H] bin. 



[Fe/H] 


Fractiono 


Factor 


Fraction 


-3.50 


0.019 


0.850 


0.016 


-3.30 


0.059 


0.790 


0.047 


-3.10 


0.140 


0.697 


0.099 


-2.90 


0.215 


0.799 


0.174 


-2.70 


0.564 


0.847 


0.484 


-2.50 


0.770 


0.845 


0.659 


-2.30 


0.642 


0.906 


0.589 


-2.10 


0.788 


0.937 


0.747 


-1.90 


1.000 


0.987 


1.000 


-1.70 


0.679 


1.005 


0.691 


-1.50 


0.782 


1.023 


0.810 


-1.30 


0.306 


0.970 


0.300 


-1.10 


0.132 


0.993 


0.133 


-0.90 


0.039 


0.999 


0.039 


-0.70 


0.205 


1.000 


0.208 


-0.50 


0.367 


1.027 


0.382 


-0.30 


0.203 


0.965 


0.199 


-0.10 


0.164 


1.022 


0.170 


0.10 


0.007 


1.069 


0.007 



To further investigate the effect of the survey-volume adjust- 
ment on the MDF, the correction procedure was also applied to 
the metal-poor giant sample of Paper V. A similar method was 
adopted, except that we assumed that the sample of Paper V are 
only giants. The corrected MDF is then compared with the ob- 
served one, as shown in the right panel of Figure [6] The sur- 
vey volume effect estimated with our method notably revises the 
shape of the giants' MDF. It clearly decreases the fraction of 
the metal-poor component and dramatically increases the pro- 
portion of the relatively metal-rich part. This effect could also 



be expected from inspection of Figure [5] because within a cer- 
tain (B - V){) bin the survey volume explored by giants with 
[Fe/H] = -3.0 (when referenced to [Fe/H] = -2.0) is obvi- 
ously larger than that of giants with [Fe/H] = -1.0, resulting 
in a much smaller correction for more metal-deficient giants. 
Thus, we conclude that although different survey volumes for 
stars with different metallicities do not affect the observed metal- 
licity distribution of a turnoff-star dominated sample, they will 
obviously change the observed MDF of a giant-dominated sam- 
ple, and cannot be ignored. In Table |H we list the correction 
factors for each [Fe/H] bin of the giants' MDF, and applied the 
values to corresponding bins of the scaled MDF of Paper V. 

3.6. Comparison with the giants' MDF 

The MDF of the HES MSTO sample can now be compared with 
that of the giants from Paper V, as shown in Figure [7] The com- 
parison between the two MDFs can be considered in two parts. 

First, at the metal-poor end with [Fe/H] < -2.0 (exclusive 
of the ultra metal-poor component with [Fe/H] < -4.0 of Paper 
V), both MDFs agree on the dramatic decrease of stars below 

[Fe/H] < -3.0 and the sharp drop at [Fe/H] 3.6. Besides, 

a ^ 2 -test of the null hypothesis that the two samples are drawn 
from the same parent distribution yields a probability of ~ 1. 
This indicates that the two samples present are quite analogous 
distributions in this metallicity region. This is not unexpected be- 
cause both the turnoff and giant samples were aimed to sample 
the Galactic halo population and were selected with similar cri- 
teria in order to derive a statistically complete sample for metal- 
poor stars. Thus, the two samples should follow similar statisti- 
cal properties in the metallicity region where the halo population 
dominates. 

The two MDFs notably differ from each other in the fraction 
of the relatively metal-rich component (e.g., [Fe/H] > -2.0), 
with the giant MDF revealing a higher fraction. For MDFs in 
the region with [Fe/H] > -2.0, the ^ 2 -test yields a probability 
of ~ 0, suggesting very different distributions. This is not diffi- 
cult to understand. As shown in the previous section, the correc- 
tion on the survey-volume has very different effects on the two 
MDFs. Also, the cutoff at (B - V) = 0.5 leads to a cutoff at com- 
paratively lower metallicites for the turnoff MDF. Hence, the two 
samples present rather distinct MDFs in this region. However, 
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Fig. 6. Comparison of the observed and the survey-volume-corrected MDF, shown with the dashed and the solid histograms, respec- 
tively. The left panel shows the HES MSTO sample, while the right panel shows the HES giant sample. 

Table 4. The correction factor arising from the survey-volume effect for corresponding [Fe/H] bins of the observed giant MDF, as 
described in Section [3~5l 



[Fe/H] 


-4.00 


-3.50 


-3.30 


-3.10 


-2.90 


-2.70 


-2.50 


-2.30 


-2.10 


-1.90 


-1.70 


Factor 


0.330 


0.331 


0.327 


0.287 


0.341 


0.407 


0.499 


0.698 


0.834 


1.184 


2.048 


[Fe/H] 


-1.50 


-1.30 


-1.10 


-0.90 


-0.70 


-0.50 


-0.30 


-0.10 


0.10 


0.30 




Factor 


3.491 


9.003 


24.97 


48.22 


107.3 


121.9 


219.0 


211.9 


2132. 


738.4 





one should keep in mind that the size of the subgroup of can- 
didates with least possibility of being metal-poor, i.e., mpcc, in 
our turnoff sample is very limited (only 16 "accepted" stars), 
and was biased against in the whole selection and observation 
procedure. Consequently, it may be incomplete for a thorough 
statistical comparison of MDFs in this [Fe/H] region. 

Therefore, as the primary motivation of this work is to dis- 
cuss the properties of the halo MDF, the completeness of both 
the turnoff and giant samples and the above quantitative inves- 
tigation should be reliable in the metallicity region which is of 
greatest interest ([Fe/H] < -2.0, especially the metal-deficient 
tail between [Fe/H] < -2.5 and -3.6). The reader should note 
that the low-metallicity tail discussed in this work is different 
from that discussed in Paper V which extends to [Fe/H] < -4.0. 

4. Do Theoretical Predictions Fit the Observations? 

One of the crucial roles that the observed halo MDF plays is to 
examine and constrain theoretical models of Galactic chemical 
evolution. In order to carry out such a comparison with any theo- 
retical predictions, we first need to convert the theoretical MDFs 
into a form that corresponds to what would be observed in a sur- 
vey with the same observational strategy and selection criteria as 
for the HES turnoff sample. 

The first modification of the theoretical MDFs is to account 
for the HES sele ction function. To a ccomplish this, we inverted 
the calibration of B eers et all (1 19991) to convert each [Fe/H] into 
a pair of (B - V)q and KP or (J - K)o and KP. Considering the fact 
that the selection function varies with (B - V)o or ( J - K)q (see 
Paper IV, V, and Section [3~4b . these theoretical "stars" were se- 
lected to follow the distribution of (B - V)o and {J - K)o of our 
observed sample. Following this, random Gaussian errors with 
standard deviations to reflect those in the measured (B - V)q and 
(J - K)q color, and the KP index were computed and added 
(c"b-v=0.06, ctj_ k =0.1, and a- KP -\.0). Finally, we applied the 
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Fig. 7. The observed MDF of HES MSTO stars (open histogram 
with solid line) is compared with that of giants by Paper V (filled 
histogram with dashed line). A relative normalization which 
scaled the maximum fraction to 1.0 is used for MDFs here and 
throughout the paper. Note that the two MDFs are scaled to 
match at [Fe/H] = -2.5 (the starting point of the metal-deficient 
tail we defined, and also marked out with the vertical dashed 
line). 



same criteria for (B — V)o or ( J - K)o versus KP to select metal- 
poor "candidates" from these theoretical stars. Using the above 
procedure we obtain a model MDF as it would have been ob- 
served in the HES (which we refer to as "as observed"). We 
compare it with the observed MDF of the turnoff sample in the 
following discussions. Since the low-metallicity tail of the MDF 
is of the greatest interest to this study, the following discussion 
will focus on the comparisons in the metallicity region between 
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[Fe/H] = -2.0 and -3.6 (where the observed MDF is considered 
statistically reliable). 

4.1. Theoretical predictions based on the Simple Model 

We begin our observational-theoretical comparison with the 
Simple Model dSearle & SargentHl972t iPagel & Patched 1 1975b 
of Galactic chemical evolution. It describes the basic form of 
a closed system which evolves from initially zero- metallicity 
as an d remains chemically homogeneous at all times. Hartwick 
1976) extended this model such that star formation ends once 
the gas is either consumed or removed (essentially relaxing the 
closure requirement of the system). Here we make use of this 
model as parameterized by the effective yield, j/ e ff, and adopting 
the same value as in Paper V, logio j/eft = -1-7. 



l.oo 



o.io 



o.oi 



Prantzos (2003) as observed 
Prantzos (2008) as observed 
HES MDP 




-2.5 -2.0 -1.5 
[Fe/H] 



1.00 



0.10 



0.01 



Hartwick (1976) as observed 
HES MDF 




-2.5 -2.0 -1.5 
[Fe/H] 



Fig. 8. The HES MS TO MDF (black histogram) is compared 
with that predicted bv lHartwickl d 1 976l) 's modified Simple Model 
(log 10 i/ e ff = -1.7) as it would have been observed in the HES 
(red solid line). Note that for all the comparisons that follow that 
we have scaled the model MDFs in such a way that they could 
best fit the metal-poor tail of the observed MDF. The vertical 
dashed line refers to [Fe/H] = -2.5 here and in the following 
plots too. 

The result is shown in Figure[8] As can be seen, the mass-loss 

modified Simple Model is able to fit the position ([Fe/H] 2.0) 

but not the height of the peak. It does, however, well fit the gen- 
eral shape of MDF tail with [Fe/H] from -2.7 through -3.6, 
although it can only predict a smooth drop of the metal-poor tail 

at [Fe/H] 3.6. This is not entirely unexpected considering the 

fact that the real Galactic halo(s) could certainly be more compli- 
cated than a simple one-zone mod el assuming th e Instantaneous 
Re cycling Approxim ation (IRA - lTinsle\fl980l) . 

IPrantzosI (l2003h addressed the effect of the IRA in the de- 
termination of the MDF of a system such as the Milky Way 
and suggested a physically motivated modification to the sim- 
ple outflow model, i.e., a composite model adopting a relaxed 
IRA, and assuming both early infall and outflow to solve the so- 
called "G dwarf problem". Based on this mod e l, and the further 
accumulation of observational data. IPrantzosI ((2008) presented 
a semi-analytical model in the framework of the hierarchical 
merging paradigm for structure formation which assumes that 
the Galactic halo is composed of the stellar debris of several sub- 
halos following either the observed properties of dwarf galaxies 
or a structure formation calculation. 



Fig. 9. Comparison of the HES M STO MD F and the models of 
Prantzos (2003) (dashed line) and lPrantzosl (120081) (solid line). 



As shown in Figure |9l both the c ompo site model with an 
early phase of gas infall by iPrant zos ( 2003) and t he hie rarchi- 
cal merging scenario for the formation by Prantzos (2008) fit the 
shape of the observed MDF tail between -2.9 and -3.4 rather 
well. However, the location of the peak of the MDF is not cor- 
rectly predicted in either case and neither of them reproduces the 
sharp drop at [Fe/H] 3.6. Rather, they predict a smooth de- 
crease of numbers of EMP stars which extend to [Fe/H] < -4.0. 



4.2. Other theoretical predictions 

Besides models based on variations of the chemical evolution 
scheme of the Simple Model, there are quite a number of other 
models based on theoretical analyses or simulations. Here we 
compare our observation with two such theoretica l predi ctions. 

The first considered is the model o f Karlsson (2006) which 
focuses on the metal-poor tail with [Fe/H] < -3.0, and attempts 
to explain the "gap" in the halo MDF with [Fe/H] between 
-4.0 and -5.0. It adoptes a scenario of negative feedback from 
Population III stars. Figure \10\ suggests that it only roughly fits 
the portion of the MDF with [Fe/H] < -3.0. It also fails to pre- 
dict the sharp drop at the low-metallicity end as well. 




-4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 
[Fe/H] 



Fig. 10. C omparison of the HES MSTO MDF and the model of 
(Karlsson 2006). 
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Another model that has been tested is GAlax y MErger Tree 
and Evolution (GAMETE - ISalvadori et"alll2007f) . It is a Monte 
Carlo code to reconstruct the merger tree of the Milky Way and 
to follow the evolution of gas and stars along the tree. This 
model defines an input parameter, the critical metallicity Z cr , 
which governs the transition from Pop III to Pop II star forma- 
tion. We compare our observed MDF with the simulated results 
corresponding to different values of Z cr , as shown in Figure [TT] 
Although according to the observational data available at that 
time, Z„. = 1O~ 4 Z was regarded as the fiducial model, it ob- 
viously cannot fit our observations here. All the predictions fail 
to fit the location of the peak of the observed MDF. Similarly to 
the conclusions in Paper V, the model with Z cr = 10~ 34 Z Q ap- 
pears to partially fit our observed MDF, being able to reproduce 
the tail with [Fe/H] < -3.0 and best predict the sharp drop at 
[Fe/H] ~ -3.6. 



1.00 



0.10 



0.01 



Z cr =10 4 a^ observed 
Z cr =10 ' ilJ as observed 
Z cr =10 ' i4 a's observed 
HES MI>£ 



-4.5 -4.0 




-2.5 -2.0 -1.5 
[Fe/H] 



0.5 



Fig. 1 1. Comp a rison of the HES MDF and the model of 
Salv adori et alj d2007l) . The purple dashed, blue dash-dotted, 
and red solid lines correspond to models with Z cr = 1(T 4 Z , 
10~ 3 5 Z Q , and 1CT 3 4 Z , respectively. 



5. Summary and Discussions 

Based on the (for now) largest metal-poor turnoff-star sample 
from the HES database and moderate-resolution follow-up ob- 
servations, we have statistically investigated the MDF of local 
MSTO stars in the Galactic halo. 

1. With reference to [Fe/H] = -2.0, the effects of relative sur- 
vey volumes have been quantitatively estimated based on 
theoretical isochrones and applied to the observed MDFs of 
both the HES turnoff and giants samples. It is shown that the 
survey-volume effect does not substantially alter the turnoff 
MDF while it dramatically changes the MDF of the giant 
sample from Paper V. 

2. The survey-volume corrected and metal-poor-class scaled 
MDFs of the turnoff sample has been compared with that 
of the halo giants. Though the two MDFs notably differ in 
the region with [Fe/H] > -2.0 (where our sample starts to 
be incomplete), for the metal-deficient region (e.g., [Fe/H] < 
-2.0) the^ 2 -test suggests that the two MDFs are quite simi- 
lar. Furthermore, both MDFs agree regarding the sharp drop 
at [Fe/H] ~ -3.6. Hence, for an MDF dominated by the halo 
population, the two MDFs agree well. 



Theoretical models of Galactic chemical evolution have been 
discussed. They can only fit portions of the observed MSTO 
MDF while none of them fully reproduces the features of the 
observations. In particular, they fail to simultaneously fit the 
peak and the metal-deficient tail between [F e/H] ~ -2.5 to 
-3.6. A lthough the Z cr = 10~ 34 Z o case of the Salvadori et al. 



(2007) model can only partially fit the observed MDF it is 

able to best predict the sharp drop at [Fe/H] 3.6. 

4. Generally, both selection criteria using KP plus (B - V)o and 
(J - K)o serve as efficient selectors of metal-poor stars. They 
are capable to reach a selection fraction up to 100% for the 
EMP candidates of our sample. 

Considering the fact that our sample mainly consists of un- 
evolved main-sequence (and subgiant) stars with low metal- 
licities, it could also provide additional useful information on 
Galactic chemical evolution. For example, a kinematic analy- 
sis of this sample could be used to re-visit the role of accretion 
of the interstellar medium during the long lifetimes of metal- 
poor stars, as approximately calcula t ed in a number of earl 
works (e.g.. iTalbot & Newmad[l977t lYoshii|[l98l Ilbenlll983l 
and also discussed by more recent studies (e.g, Christlie b et al 
120041: iNorris et al.ll2007tlFrebel et al.ll2009l) . 

It should also be pointed out that all of our comparisons of 
the MDFs have been performed under the assumption that we 
are modeling a single halo population, which current evidence 
suggests is an over-simplification. It seems likely that the ob- 
served MDFs for both the HES MSTO stars and the HES giants 
comprise overlapping contributions from the outer-halo popu- 
lation at the lowest metallicities and the inner halo at interme- 
diate low metallicities, with respective tails of as-yet unknown 
relative strengths and convolved with the HES metallicity se- 
lection bias that becomes more severe above [Fe/H] ~ -2.5 to 
-2.0. This possibility was already mentioned in Paper V where 
it was noted that there appeared to be relatively larger fractions 
of EMP stars at heights above the plane |Z| > 15 kpc than in the 
intermediate range 5 < |Z| < 15 kpc. This is in line with t he ex- 
pecta tions of the dual halo interpretation of ICarollo et alJ ({2007, 
2010). Progress on this issue will come from consideration of 
the dual halo modeling approach, ideally in combination with a 
full kinematic analysis of these samples that forms the basis of a 
paper in preparation. 

However, the HES metal-poor turnoff sample discussed in 
this paper contains no objects with [Fe/H] < -3.6 which obvi- 
ously do exist. Thus, we are not able to discuss the performance 
of theoretical MDFs in the most metal-deficient regime. Larger 
statistically complete samples are required for a thorough com- 
parison with theoretical predictions. Fortunately, such samples 
will be obtained from much larger and deeper surveys in the near 
future, such as from SEGUE-2 and the Apache POint Galactic 
Evolution Experiment (APOGEE), the Large Sky Area Multi- 
Objec t Fiber Spectroscopic Telescop e (LAMOS T, Zhao et all 
l2006h . and the Southern Sky Survey (keller et al.ll2007l) . 
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